PCAC and coherent pion production by low energy neutrinos 



Ch. Berger * 

I. Physikalisches Institut dcr RWTH Aachen University, Germany 

L. M. Sehgal ^ 

Institut fiir Theoretische Physik (E) der RWTH Aachen University, Germany 



Abstract 

Coherent tt"*" and 7r° production in low energy neutrino reactions is discussed in the 
framework of the partially conserved axial vector current theory (PCAC). The role of 
lepton mass effects in suppressing the 77+ production is emphasized. Instead of using 
models of pion nucleus scattering, the available data on pion Carbon scattering are 
implemented for an analysis of the PCAC prediction. Our results agree well with the 
published upper limits for tt"*" production but are much below the recent MiniBooNE 
result for 7r° production. 

PACS 13.15.+g 11.40.Ha 25.80.Dj 

1 Introduction 

The availability of high intensity neutrino beams with energies up to a few GeV opens 
the way to precise investigation of neutrino oscillations. Essential for these experiments 
is a detailed understanding of all low energy neutrino reactions especially single pion 
production in charged (CC) and neutral (NC) current reactions. Coherent pion production 
off nuclei, e.g. in C — > z^^ C + vr" constitutes an especially interesting subsample 

not only because it is a significant background to the — > z^e oscillation search but also 
because it is deeply rooted in fundamental physics via Adler's PCAC theorem [HE] which 
connects forward neutrino scattering with the pion nucleon cross section. 

2 PCAC and forward lepton theorem 

Our starting point is the general expression for inelastic neutrino scattering 

WdTy = 4vr^ """^ [ua, + van + 2uvas] (1) 
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^Gf and 6c are the Fermi coupling constant and the Cabbibo angle. 
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already derived by Lee and Yang j3] in 1962 for zero mass of the outgoing lepton. The 
momentum and energy transfer between incoming neutrino and outgoing lepton is given 
by q and u = E — E' with y = u/E. As usual = —q^ denotes the four-momentum 
transfer squared, = — while k = (W"^ — Mfj)/2Mj\[ for a hadronic system with 
invariant mass W emerging from a nucleus or nucleon with mass M^. The kinematical 
factors u, v are given hy u,v = (E + E' i: \q\)/2E. 

The cross sections crL,R,s (where L, R, S stands for Left, Right, Scalar) are unknown 
functions of and W which have to be measured or calculated in theoretical models. It 
is well known that for ^ only the term with as in ([T|) survives because the overall 
factor is compensated by a factor 1/Q^ in the scalar cross section. In this limit PCAC 
predicts [1] 

Here /tt is the pion decay constant (130.7 MeV) and (T,rAr(VF) the pion nucleon (or nucleus) 
cross section for the hadronic final state under consideration. 

The resulting formula for forward inelastic neutrino scattering i^ 



dQ^dy 



G'p cos' 9c fi E 



uva^NiW) . (3) 



In the presence of lepton mass nii there is a correction to this formula caused by the pion 
pole term in the hadronic axial vector current. Including, in addition, an axial vector form 
factor Ga{Q'^) to describe the variation of the cross section at small values of around 
the forward direction, one obtains the result given by Kopeliovich and Marage [5] 



da^^ Gl cos' Ocf^E 

— -- — rUV 



dQ'dy 27r2 



ri ^ Qmin A I y (n2 \ Qmin 



2Q' + mlJ 4'^ ^'""^'(g2 + ^2)2 

(4) 

where Q^j^ = mfy/{l — y) is the high energy approximation to the true minimal Q^. The 
axial vector form factor Ga is defined by 

with a typical value for the axial vector meson mass ttia of 0.95 GeV; see however the 
extensive discussion in the literature [7]. 

The first term inside the rectangular brackets of @ corresponds to outgoing muons 
with negative helicity (helicity nonflip) whereas the second term is the helicity flip con- 
tribution which vanishes at 0° scattering angle. This helicity structure was also found 
by Adler [2] and Piketty and Stodolsky [6]. With Ga = 1 the expression inside the rect- 
angular brackets represents the Adler screening factor which was invoked in [9l [10] as a 
possible explanation of the dip in CC reactions at low Q', resulting from the destructive 
interference of the pion pole. 



Oi denotes the laboratory angle of the outgoing lepton. 
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For neutrino scattering off nuclei N the coherent pion channel f^iV — > h^tt'^N has 
special interest, since the PC AC formula ^ predicts the cross section to be proportional 
to the elastic cross section vrA^ — > vrA^. This hadronic process is strongly enhanced in the 
forward direction because of the coherent action of the A nucleons in the nucleus, with an 
amplitude ~ A at 0° scattering angle. This effect then implies the enhancement of forward 
going pions in the neutrino reaction, which is a characteristic signature of coherent pion 
production, distinguishing it from possible incoherent backgrounds. An early analysis of 
pion production by neutrinos, that included a discussion of the kinematical region termed 
coherent, was given in [8]. 



3 Application to coherent pion production 



In the following we investigate in more detail coherent single pion production. Assuming 
that the derivation given above also holds for the differential cross section one gets for the 



CC reaction u^N 



dQ^dydt 



27r2 



uv 



)2. 

' mm 



2 Q2 + ^2 



Q 



ram) 



Q 



2 

min 



(g2 + ^2)2 



da{Tr+N tt+N) 



dt 



(6) 



where t is the modulus of the four-momentum transfer squared between incoming virtual 
boson and outgoing pion. For calculating the NC reaction vN — > vk°N one has to set 
= 0, 6c = Q and divide the right hand side of the resulting equation by 2, because 
f-K° = /tt/ \/2. We thus obtain 



da 



NC 



■uvG 



dQ'^dydt 4^2 |g 



dt 



(7) 



For isoscalar targets da{iT^N — > vr+A^) equals da{'K°N — > 7r°N). 

In the widely used Rein-Sehgal (RS) model [TT] the kinematical factors on the right 
hand side of ([7]) have been evaluated for = resulting in the simpler expression 



da^^ _Glfll 



dQ^dydt 4tt^ y 



/ 2 da{T:°N 

-Li A 



7r°A^) 



dt 



(8) 



for the cross section. 

We have compared ([7]) to ([8]) by performing a Monte Carlo integration using the 
illustrative ansatz 



dt 



aobe 



(9) 



with energy independent coefficients ctq = 80 mb and & = 45 GeV"^ for elastic 7r° or vr"*" 
scattering off Carbon nuclei. (This ansatz will be motivated in the next section.) 

The result of the integration is shown in figH^ where the ratio Cfjju/cr^jjjjpig is plotted 

versus the energy of the incoming neutrino. Here fi^^ stands for the integral of ([7]) and 
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Figure 1: a) Ratio o■J[^^/c^^i|!f^pl(, of the integrated cross sections of ([7]) and ([5]) versus the 
energy of the incoming neutrino, b) Ratio a'^^ /a'"" of the integrated cross 
sections (O and ((T]) versus the energy of the incoming neutrino. 



"^simple ^^"^ integral of This figure demonstrates that at low energies the neutrino 
cross section is substantially reduced by using the complete kinematical factors of 

We also compared tt"*" production to 7r° production by integrating ([6]) and ([7]) using 
identical form factors and pion nucleus cross sections. Apart from the factor cos^ Oq one 
would naively expect a ratio of 2 according to the ratio of the pion decay constants. FigHb 
shows however a remarkable violation of this isospin symmetry due to lepton mass effects 
contained in ([6|). The dominant cause for the variation observed in the figure is the reduced 
phase space of CC reactions. The Adler screening factor reduces the CC cross section by 
further 10% at ii^ = 0.6 GeV and 4% at E = 2 GeV. The influence of the numerical value 
of niA on the cross section ratio is neglegible. Even setting Ga = 1 changes the ratio by 
less then 2%. 

4 The elastic pion Carbon cross section 

In [11] a model has been presented which calculates elastic pion nucleus scattering from 
pion nucleon scattering via 

da{7rN ^ ttN) ^2'^'^ci 
di ~ "dT 

Here the pion on the left hand side can be charged or neutral. The differential elastic pion 
nucleon cross section in forward direction on the right hand side is determined via the 
optical theorem (neglecting a possible real part of the scattering amplitude) 



-bt : 



abs 



(10) 



and the slope b of the exponential t-distribution is taken from the optical model relation 

b = ^RlA'/' (12) 

with e.g. Rq = 1.057 fm. Fabs describes the average attenuation of a pion emerging from 
a sphere of nuclear matter with radius RqA^^^ resulting in 
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Fabs = exp I - -^Q^j^2 ^incl ) (^^) 



with 



^inel = • (14) 



As an example we calculate the total elastic pion Carbon cross section via 



aei(7r C ^ TT C) = I I . (15) 



The total pion nucleon cross sections are available as computer readable files [12]. The 
data were fitted by a superposition of Breit Wigner functions and a Regge inspired term 
+ oi/y^Ip^I with Ip^I denoting the pion laboratory momentum. A similar fit to the 
elastic cross sections finally yields (Jinei = Ctot — which is used for calculation of Fabs- 

The dotted line of figl2] shows the result. The fact that this cross section is derived from 
a simple (classical) ansatz expressed by (fTOj) and (fT3|) raises doubts about its validity as a 
description of pion-nucleus scattering in the resonance region. An alternative approach to 
coherent pion nucleus interaction based on the Glauber model was proposed by Bel'kov 
and Kopeliovich [13]. Its numerical results were similar to those in the RS model at least 
at high energies. The experimental groups use detailed Monte Carlo routines to simulate 
the scattering and absorption of the pion inside the nucleus |14t [T^ . 

Following the PCAC route we have tried to circumvent the uncertainties in modelling 
nuclear processes by direct appeal to data on pion nucleus elastic scattering, see also [TB] . 
For Carbon targets this can be done easily because vr"*" and vr" data on differential and 
total cross sections exist for pion kinetic energies from 30 to 870 MeV. They have been 
subjected to phase shift analyses yielding up to 21 complex phase shifts per energy [T7] . 
Neglecting electromagnetic effects these phase shifts can be used to compute the elastic 
strong interaction cross section dcTd/dt for pion Carbon scattering in a straightforward 
manner. 

The phase shifts accurately reproduce even tiny effects like secondary peaks in the 
angular distribution. We have checked that except for the lowest two kinetic energies of 
30 and 50 MeV it suffices to parametrize the cross section by the simple ansatz 



dt 



Aie""^* (16) 



with energy dependent coefficients , 6i , which are listed in the table on the right hand 
side of fig|2j For energies between the measured data points these coefficients are linearly 
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Figure 2: Total elastic pion Carbon cross section versus pion laboratory momentum. 

The dotted line represents the Rein-Sehgal model according to (fT5|) . the sohd 
line is derived from pion Carbon data as explained in the text. The table on 
the right hand side contains the coefficients Ai, hi of (|16p . 



interpolated which is the reason for the zig-zag structure of the solid line in fig|2j It is 
obvious that (TcI from pion Carbon data is much below the RS model in the resonance 
region. At the same time one observes that as \p^\ approaches 1 GeV, the two curves 
become very similar with a^x ~ 80 mb. This finally justifies the ansatz Q. It also 
suggests that the RS hadronic model fails in the region of the A resonance, but may be a 
valid description at higher energies. 

5 Results 

We are now ready to integrate the cross section ([6]) for the two different models of pion 
Carbon scattering discussed in the last section. The results are plotted versus the neutrino 
energy in figlSti for vr'^ production and in figlHb for vr"^ production. In obtaining the 
lower curves the empirical pion Carbon cross sections were calculated by assuming the 
coefficients in the last line of the table on the right hand side of figi2] to be valid up to 
= 1.7 GeV. An error of 30% in this assumption results in a cross section error of 6% 
s± E = 2 GeV. 

The curve using Carbon data is a factor of 3 to 2 below the curve obtained by applying 
the RS hadronic model. Cross sections for NC and CC coherent single pion production 
on Carbon have also been calculated using an ansatz based mainly on the microscopic 
process vp /i^A'^"*" and its modification in the nuclear environment [181 [T9l [20| 121] . 
(For an early reference to this subject see [22].) Remarkably our calculations agree well 
with the corresponding results given in |20[ |2I] based on a very different approach to 
coherent neutrino scattering. The predicted cross sections of [l6] depend sensitively on a 
cut paarameter ^. Referring to footnote 41 of [TB] with ^ = 1 the results are close to the 
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Figure 3: Cross section per nucleus of colierent tt production by neutrinos off Carbon 
nuclei, a) NC reaction i/^ ^ _^12 ^ qq reaction C -> 

/x~ C + 7r+. The results in units of 10"'"' cm^ are plotted versus the 
neutrino energy in GeV. The upper curve is calculated using the hadronic RS 
model, the lower curve using our parametrization of pion Carbon scattering 
data. 




Figure 4: Differential cross section da/dQ^ per nucleus for coherent single pion produc- 
tion off Carbon nuclei. The data are obtained by integrating ([7]) using Carbon 
data for cr,rjv- The neutrino energy is 1 GeV. Solid line is for the NC reaction, 
the dashed line for the CC reaction. 
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ones obtained in this paper. The differential cross sections da/dQ"^ or da/dcos 9i are more 
sensitive to details of the theoretical models. We give in figUlour prediction for da/dQ"^ 
at a neutrino energy of 1 GeV. For the CC reaction a pronounced dip in forward direction 
is seen which is mainly due to the Adler screening factor contained in the rectangular 
brackets of ([7]). 

Comparing to experimental results we first discuss vr"*" production. At a neutrino 
energy of 1.3 GeV the K2K experiment [23] has placed an upper limit of 0.60 x 10~^ on 
the cross section ratio of coherent pion production cr^^^ to the total CC cross section cr^*^ 
at 90% confidence level (CL). Using the K2K cut in the muon momentum (p^ > 0.45 GeV) 
our prediction is 0"^^^ = 0.62 x 10"^'' cm^ per nucleon. With a^^ = 107 x 10~^^cm^ as 
quoted in [23j we obtain cr^oh/'^?'^ ~ ^-^^ ^ 10~^ which is consistent with the K2K result. 

The SciBooNE experiment [K] has measured an upper limit of 0.67 x 10~^ at 1.1 GeV 
(90% CL) on the same cross section ratio. Assuming the a'^^ value quoted in [23] to scale 
linearly with E we obtain a'^^ /E = 82.3 x lO""^*^ cm^/GeV in the low energy regime and 
al^JcTu^ = 0.58 X 10-2 taking al^^ from figEb. At 2.2 GeV the SciBooNE upper hmit 
is 1.36 X 10"^ to be compared with a calculated ratio of 0.68 x 10"^. The latter value 
requires a slight extrapolation of figlSb- 

The improved model for coherent pion production is thus in good agreement with the 
TT^ measurements. It is possible that the lower cross section for coherent tt"'' produc- 
tion also reduces the discrepancy at small in CC reactions noted by the MiniBooNE 
collaboration |24j . 

Turning to 7r° production the published coherent fraction (Tcoh/(c"coh + Cincoh) of the 
MiniBooNE experiment [H] is (19.5 it l.l(stat) ib2.5(sys))% at a nominal neutrino energy 
of 1.2 GeV. Using the Rein-Sehgal model for incoherent production [25] and our present 
calculation of the coherent cross section we get a coherent fraction of « 5% which is much 
below the experimental findings. It would be interesting to see if there is a dependence 
of the experimental fraction on the details of the coherent model used in the analysis. 
It should be stressed that the theoretical prediction for coherent scattering covers only 
reactions where the nucleus stays intact and does not break up during interaction. The 
data on pion-Carbon scattering [T7] show that nearly two-thirds of the cross section is 
inelastic. Our estimate of the error in the theoretical prediction (taking account of the 
model dependence of Uincoh; the extrapolation in Q^, the interpolation in fig|2] and the 
neglect of the transverse contributions (Tr^l estimated in [16]) is 20%. 

The total pion nucleus cross section scales ~ A^^^ at least for small atomic numbers 
A |26j . If the ratio of elastic to total cross section depends only weakly on A we then 
expect the elastic cross section also to scale proportional to A^/^ . Using the RS hadronic 
model the ratio of the cross sections for Aluminum [A = 27) and Carbon is 1.68 for E = 2 
GeV corresponding to a scaling law ~ A^'^^ . The algorithm presented in the preceding 
section can thus probably be simply extended to other light nuclei. Applying the A^^^ 
scaling law we obtain for Aluminum cr^^j^ = 13.3 x lO"^'^ cm^ at £" = 2 GeV which agrees 
within errors with the experimental value of (27zb 17) x lO"^'^ cm^ measured by the Aachen 
Padua experiment [27] . 

For neutrino energies > 2 GeV the calculation using Carbon data is not applicable 
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because a large part of the integration then needs Carbon data for \p^\ > 1 GeV which 
are not available. The fact that the RS prediction for pion Carbon scattering overlaps 
with the empirical result in the region \p^\ > 0.7 GeV (figl2|) may be a hint that the RS 
model begins to be a valid description of coherent pion production beyond the resonance 
region. Such a view would explain the impressive agreement of the RS model with the 
large body of data on coherent pion production in the energy domain E = 2. .100 GeV, as 
documented for example in [28]. 
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